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An analysis of inclusive quasielastic electron scattering is presented using different descriptions 
of the final state interactions within the framework of the relativistic impulse approximation. The 
relativistic Green's function approach is compared with calculations based on the use of relativistic 
purely real mean field potentials in the final state. Both approaches lead to a redistribution of the 
strength but conserving the total fiux. Results for the differential cross section at different energies 
are presented. Scaling properties are also analyzed and discussed. 
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I. INTRODUCTION 

Electron scattering reactions with nuclei have provided 
the most detailed and complete information on nuclear 
and nucleon structure. Analysis of data for light-to-heavy 
nuclei and for different kinematical situations have been 
presented in the literature [1-5]. Not only differential 
cross sections but also the contribution of the separate 
response functions have been considered. From the the- 
oretical point of view, an important effort has been de- 
voted for the last 20 years to the description of inclu- 
sive and exclusive processes. The high energies and mo- 
menta involved in recent experiments have led to the 
development of fully relativistic models describing the 
scattering process [6-17]. Moreover, ingredients beyond 
the quasielastic (QE) approach, namely, meson exchange 
currents, correlations, etc., have been also considered fol- 
lowing different approaches [18-32]. 

Within the QE domain, which is the region consid- 
ered in this work, the treatment of final state interactions 
(FSI) between the ejected nucleon and the residual nu- 
cleus has been proved to be essential in order to compare 
with data. This has been well established in the case 
of exclusive (e, e' N) reactions where an ejected nucleon 
is detected in coincidence with the scattered electron. 
Analyses based on the use of phenomcnological complex 
potentials in the final channel have provided results in 
accordance with data [1, 4]. In particular, the use of 
relativistic complex optical potentials within a fully rel- 
ativistic description of the reaction mechanism has led 
to theoretical cross sections in excellent accordance with 
data [7, 8, 10, 11, 13, 15, 16]. Moreover, comparison with 
separate responses and asymmetries has also proved the 
capability of the relativistic approach to successfully de- 
scribe fine details of data behavior [9, 10, 13, 33, 34]. 

In the analysis of inclusive reactions, contrary to exclu- 
sive ones, all inelastic channels in the final state should 
be retained. This means that the flux is conserved, and 



consequently, the distorted wave impulse approximation 
(DWIA) based on the use of a complex potential should 
be dismissed. However, final state interaction continues 
to be a main ingredient in the inclusive process, and its 
appropriate description is required in order to describe 
data. Within the framework of fully relativistic models 
different approaches have been presented in the litera- 
ture. On the one hand, a description based on the rela- 
tivistic distorted wave impulse approximation (RDWIA) 
has been pursued, but using purely real potentials in the 
final channel. This is consistent with flux conservation. 
Concerning the specific potentials, the final nucleon state 
has been evaluated with the real part of the relativistic 
energy-dependent optical potential, denoted as rROP, or 
with the same relativistic mean field potential consid- 
ered in describing the initial nucleon state, denoted as 
RMF (see [15-17, 35]). A second approach of FSI makes 
use of the relativistic Green's function technique, where 
the components of the nuclear response are written in 
terms of the single particle optical model Green's func- 
tion. This method, denoted as GF, allows one to perform 
calculations treating FSI consistently in the inclusive and 
exclusive channels. The same relativistic (complex) op- 
tical potential is considered in both cases, but flux is 
conserved in the inclusive process. Moreover, redistribu- 
tion of strength among different channels is due to the 
real and also signiflcantly to the imaginary part of the 
potential (see [14, 36-39] for details). 

The exhaustive analysis of the (e, e') world data has 
demonstrated the quality of scaling arguments at high 
momentum transfer for excitation energies below the QE 
peak [40-43] . These data, and particularly those coming 
from the separate longitudinal contribution, have been 
shown to respect scaling of first kind (independence of 
the transferred momentum q) and scaling of second kind 
(no dependence on the nuclear system). The analysis of 
the longitudinal response, very mildly affected by me- 
son exchange currents and nuclear correlations, has per- 
mitted the extraction of a phenomcnological QE super- 
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scaling function given for transferred energies below and 
above the QE peak. These regions correspond to neg- 
ative and positive values of the so-called supcrscaling 
variable, respectively [35, 41-44]. Scaling analysis for 
electron scattering, which was extended into the A re- 
gion [20, 45, 46], has been also exploited to predict in- 
clusive charged-current (CC) neutrino- nucleus cross sec- 
tions [35, 44, 45, 47-50] and neutral-current (NC) pro- 
cesses [51-53]. 

The QE phenomenological scaling function presents a 
significant asymmetry with a tail extended to high val- 
ues of the transferred energy u> (positive values of the 
scaling variable tjj'{q,u))). This imposes strong restric- 
tions to all theoretical models describing QE (e, e') pro- 
cesses, namely, they should be able to fulfill scaling prop- 
erties and in addition reproduce the specific shape of 
the scaling function. Asymmetry in the scaling func- 
tion is largely absent in non-relativistic models based on 
a mean field approach [48]. The same comment applies 
to the rclativistic plane wave limit and to results based 
on a DWIA using real rclativistic energy-dependent op- 
tical potentials [35, 44]. In contrast, asymmetry comes 
out within the rclativistic impulse approximation, but 
with FSI described using strong rclativistic mean field 
potentials. An asymmetrical scaling function has been 
also shown to occur within the framework of a semi- 
relativistic model, based on FSI given through the Dirac 
equation-based (DEB) potential derived from the RMF 
(see [54] for details). 

Previous studies clearly illustrate the central role 
played by FSI in providing theoretical results in accor- 
dance with data. Hence, in this work we present a care- 
ful analysis of (e,e') processes comparing two basic de- 
scriptions of FSI: the RMF approach and the rclativistic 
GF technique. Both models have their own merits. In 
the RMF the flux lost into inelasticities is recovered by 
using the same, purely real, energy-independent poten- 
tial seen by the bound nucleons and thus no information 
from scattering reactions is explicitly incorporated in the 
model. In this way it is consistent with dispersion rela- 
tions [55]. The RMF leads to an asymmetrical scaling 
function which is supported on overall by data behavior. 
However, it breaks down scaling significantly as the trans- 
ferred momentum q increases, particularly in the region 
above the QE peak. The analysis of experimental (e, e') 
data indeed leaves room for first-kind scaling breaking in 
this region, due partly to A production and other contri- 
butions beyond the impulse approximation, namely, me- 
son exchange currents and their impact in the 2p—2h sec- 
tor [24]. The GF approach provides a consistent and uni- 
fied treatment of inclusive and exclusive electron-nucleus 
scattering processes. It also fuUfiUs dispersion relations 
and recovers the flux lost into inelasticities by means of a 
formal summation of unobserved channels performed via 
the optical potential. Most often, the optical potential is 
taken from phenomenology, and thus, in this approach, 
information on the observed nucleon-nucleus scattering 
is incorporated into the model. 



The differences observed between the predictions of the 
two approaches and their behaviour with regard to scal- 
ing arguments will be helpful in settling down an ap- 
propriate description of FSI for inclusive (e, e') processes 
at different kinematics. Further, this would be useful for 
understanding nuclear effects, especially FSI, in neutrino- 
nucleus cross-sections [35, 37, 44, 45, 56-65], which are of 
paramount interest for the ongoing and future neutrino 
oscillation experiments [66 74]. 

Therefore, in this work we present a systematic study 
of both models, comparing their predictions and analyz- 
ing their scaling properties, with special emphasis on the 
specific shape of the scaling function. 

The paper is organized as follows. In Sec. II, 
we present the basic formalism involved in inclusive 
quasielastic electron scattering with a discussion of the 
different models considered in the description of FSI, 
namely, the rclativistic mean field and the rclativistic 
Green's function approach. We also provide a subsection 
which contains the basic facts related to scaling argu- 
ments and introduce supcrscaling functions. In Sec. HI, 
we show and discuss our results. First, we focus on the 
analysis of differential cross sections and later on we an- 
alyze in detail results for the scaling function. Finally, in 
Sec. IV we summarize the main results and present our 
conclusions. 



II. INCLUSIVE QUASIELASTIC ELECTRON 
SCATTERING 



An electron with four-momentum /ST^ = (e, k) is scat- 
tered through an angle ??e to four-momentum K'^ = 
{e',k'). The four-momentum transfer is = — 
K'^ = {LO,q). In the one-photon exchange approxima- 
tion the inclusive cross section for the quasielastic (e, e') 
scattering on a nucleus is given by [1] 



ds' dn' 



^aMlVLRL + VTRr] , 



(1) 



where f2' is the scattered electron solid angle and ctm is 
the Mott cross section [1]. The coefficients V come from 
the lepton tensor components and are defined as 



Vrp = tan' 



2|g|' 



(2) 



where \Q'^\ = jqp -w^. 

All nuclear structure information is contained in the 
longitudinal and transverse response functions, Rl and 
Rt, expressed by 



RL{q,uj) 
RT{q,oj) 



W^\q,uj) , 



(3) 



in terms of the diagonal components of the hadron tensor 

that is given by bilinear products of the transition ma- 
trix elements of the nuclear electromagnetic many-body 
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current operator between the initial state | ^fo) of the 
nucleus, of energy Eq, and the final states | ^l/f), of energy 
Ei, both eigenstates of the (A+ l)-body Hamiltonian H, 
as 

X (*o I rHq) I */) m^ + uj-Ef),{A) 

involving an average over the initial states and a sum over 
the undetected final states. The sum runs over the scat- 
tering states corresponding to all of the allowed asymp- 
totic configurations and includes possible discrete states. 

In the QE region the nuclear response is dominated 
by one-nucleon knockout processes, where the scatter- 
ing occurs with only one nucleon which is subsequently 
emitted. The remaining nucleons of the target behave 
as simple spectators. Therefore, QE electron scattering 
is adequately described in the Relativistic Impulse Ap- 
proximation (RIA) by the sum of incoherent processes 
involving only one nucleon scattering. 

In the RIA the nuclear current operator is assumed to 
be the sum of single- nucleon currents j^^ , for which differ- 
ent relativistic free nucleon expressions [75] can be used. 
In the present calculations we use the option denoted as 
CC2, i.e., 

f = F,{Q^)r+i^F2{Q^)a^''Q., (5) 

where k is the anomalous part of the magnetic moment 
and (7^"' = {i/2) [j'',^'']. Fi and F2 are the Dirac and 
Pauli nucleon form factors [76]. 

Within the RIA framework and under the assumption 
of a shell-model description for nuclear structure, the 
components of the hadron tensor are obtained from the 
sum, over all the single-particle (s.p.) shell-model states, 
of the squared absolute value of the transition matrix 
elements of the single-nucleon current 

{x^^-\E)\riq)\^^) , (6) 

where Xe'\^) the scattering state of the emitted nu- 
cleon and the overlap ipn between the ground state of 
the target | ^fo) and the final state | n) of the residual 
nucleus is a s.p. shell-model state. 

In the calculations presented in this work the bound 
nucleon states arc taken as self-consistent Dirac- 
Hartree solutions derived within a RMF approach using 
a Lagrangian containing a, oj and p mesons [77-81]. 

Different prescriptions are used to calculate the rel- 
ativistic s.p. scattering wave functions. In the simplest 
approach the plane- wave limit is considered, i.e., the Rel- 
ativistic Plane- Wave Impulse Approximation (RPWIA), 
where FSI between the outgoing nucleon and the resid- 
ual nucleus arc neglected. In the approaches based on the 
relativistic distorted- wave impulse approximation FSI ef- 
fects are accounted for by solving the Dirac equation 
with relativistic optical potentials. The use of relativis- 
tic energy-dependent complex optical potentials fitted to 



elastic proton-nucleus scattering data has been very suc- 
cessfuU in describing exclusive (e,e'p) data [1, 4, 7 13]. 
In the exclusive scattering the imaginary part of the op- 
tical potential produces an absorption that reduces the 
calculated cross section and accounts for the flux lost 
in the considered elastic channel towards other channels. 
This approach is appropriate for the exclusive process 
where only one channel is considered, but it would be 
inconsistent for the inclusive scattering, where all the in- 
elastic channels should be retained and the total flux, 
although redistributed among all possible channels due 
to FSI, must be conserved. As a result, in the RDWIA 
(with complex potentials) the flux is not conserved and 
the inclusive (e, e') cross section is underestimated [15- 
17, 36]. A simple way of estimating the inclusive re- 
sponse within the RIA, avoiding spurious flux absorp- 
tion, is to use purely real potentials. In a first approach, 
the imaginary part of the phenomenological relativistic 
energy-dependent optical potentials [82] is set to zero, 
thus retaining in the calculations only the real part. In 
a second approach, the scattering states are described 
by using the same RMF theory potentials considered in 
the description of the initial bound state (p„. We refer 
to these two different FSI descriptions as real Relativis- 
tic Optical Potential (rROP) and Relativistic Mean Field 
(RMF), respectively. 

We note that rROP conserves the flux and thus it is 
inconsistent with the exclusive process, where a complex 
optical potential must be used. Moreover, the use of a 
real optical potential is unsatisfactory from a theoreti- 
cal point of view, since it is an energy-dependent poten- 
tial, reflecting the different contribution of open inelastic 
channels for each energy. Dispersion relation then dic- 
tates that the potential should have a nonzero imaginary 
term [55]. On the other hand, the RMF model is based 
on the use of the same strong energy-independent real 
potential for both bound and scattering states, so that it 
fulfills the dispersion relation [55] and, further, it main- 
tains the continuity equation. 

Green's function techniques are exploited to derive the 
inclusive response in a different model where the flux is 
conserved and the use of a complex optical potential al- 
lows one to treat FSI consistently in the inclusive and 
in the exclusive reactions. Detailed discussions of this 
Green's function (GF) approach can be found in [14, 36 
39]. Here we recall only the essential features of the 
model. 

In the GF approach the components of the nuclear re- 
sponse in Eq. (4) are written in terms of the s.p. optical 
model Green's function. This is the result of suitable 
approximations, such as the one-body current assump- 
tion and subtler approximations related to the lA. The 
explicit calculation of the s.p. Green's function can be 
avoided by its spectral representation, which is based on a 
biorthogonal expansion in terms of a non Hcrmitian opti- 
cal potential and of its Hermitian conjugate. The nuclear 
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response of Eq. (4) is then obtained in the form [36] 



ReTi:''{Et-en,Et-en) 



--V r d£- i -Imrr(f,£;f-en) 



(7) 



where n denotes the cigcnstate of the residual nucleus 
related to the discrete eigenvalue £„ and 



Tit<^{e,E) = {^n\f\q)yJl-V{E)\x^i\E)) 

X (x^-)(£) I sJl-V{E)f{ci) I Vn) .(8) 



The factor ^1 - V{E), where V {E) is the energy 
derivative of the optical potential, accounts for interfer- 
ence effects between different channels and justifies the 
replacement in the calculations of the Feshbach optical 
potential V of the GF model, for which neither micro- 
scopic nor empirical calculations are available, by the 
local phenomenological optical potential [36, 38]. Dis- 
regarding the square root correction, the second matrix 
element in Eq. (8) is the transition amplitude of single- 
nucleon knockout of Eq. (6), where the imaginary part 
of the optical potential accounts for the flux lost in the 
channel n towards the channels different from n. In the 
inclusive response this loss is compensated by a corre- 
sponding gain of flux due to the flux lost, towards the 
channel n, by the other final states asymptotically origi- 
nated by the channels different from n. This compensa- 
tion is performed by the first matrix element in the right 
hand side of Eq. (8), that is similar to the matrix ele- 
ment of Eq. (6) but involves the eigenfunction x^^\e) 
of the Hermitian conjugate optical potential, where the 
imaginary part has an opposite sign and has the effect of 
increasing the strength. In the GF approach the imag- 
inary part of the optical potential redistributes the flux 
lost in a channel in the other channels, and in the sum 
over n the total flux is conserved. 

The hadron tensor in Eq. (7) is the sum of two terms. 
The calculation of the second term requires integration 
over all the cigenfunctions of the continuum spectrum of 
the optical potential. If the imaginary part of the optical 
potential is neglected, the second term in Eq. (7) vanishes 
and; but for the square root factor (whose contribTition 
is however generally small in the calculations), the first 
term gives the rROP approach. 



A. Scaling at the queisielastic peak 

Scaling ideas applied to inclusive QE electron-nucleus 
scattering reactions have been shown to work properly 
to high accuracy [41-43]. In fact, an exhaustive analysis 
of QE (e, e') world data has demonstrated that scaling 
property is well respected at high momentum transfer 
for excitation energies falling below the QE peak. These 



data, when plotted against a properly chosen variable 
(the scaling variable), show a mild dependence on the 
momentum transfer (reasonable scaling of first kind) and 
almost no dependence on the nuclear target (excellent 
scaling of second kind). Simultaneous fulfillment of both 
conditions leads to superscaling. The analysis of the sep- 
arated longitudinal contribution to (e,e') data has led 
to introduce a phenomenological scaling function which 
contains the relevant information about the nuclear dy- 
namics explored by the probe. Scaling and superscaling 
are general phenomena exhibited by nature which any 
"reliable" model providing a description of the scattering 
process should be able to reproduce. Not only superscal- 
ing behavior should be described, but also the specific 
magnitude and shape of the "universal" phenomenologi- 
cal superscaling function needs to be reproduced. 

The usual procedure to get the scaling function con- 
sists of dividing the inclusive differential cross section (1) 
by the appropriate single-nucleon eN elastic cross sec- 
tion weighted by the corresponding proton and neutron 
numbers [41-43, 83] involved in the process, 



da 



de'dn' 



(TM [VLGL{q,u)) + VTGT{q,u))] ' 



(9) 



Analogously, the analysis of the separated longitudinal 
(L) and transverse (T) contributions leads to scaling 
functions. 



THT){->P ,q) = kp-^ J r 



(10) 



The term ?^'(g,a;) is the dimensionless scaling variable 
extracted from the relativistic Fermi gas (RFG) analysis 
that incorporates the typical momentum scale for the se- 
lected nucleus [35, 41]. The fully relativistic expressions 
for Gi(T) involve the proton and neutron form factors 
'-^B(M)' weighted by the proton and neutron numbers, 
and an additional dependence on the nuclear scale given 
through the Fermi momentum kp (explicit expressions 
are given by Eqs. (16-19) in [41]). 

Whereas L data are compatible with superscaling be- 
havior, permitting the extraction of the phenomenolog- 
ical function f^^{fp'), scaling is known to be violated 
in the T channel at energies above the QE peak by ef- 
fects beyond the impulse approximation [20 24, 84]. It 
is important to point out that, although many models 
based on the impulse approximation exhibit superscal- 
ing, even perfectly as the RFG, only a few of them are 
capable to accurately reproduce the asymmetric shape of 
f^P with a significant tail extended to high transferred 
energies (large positive values of the scaling variable tp'). 
One of these is based on the RIA with final state interac- 
tions described by means of strong relativistic mean field 
potentials. On the contrary, calculations based on the 
plane wave limit and/or the use of real relativistic optical 
potentials in the final state, although satisfying super- 
scaling properties, lead to "symmetrical-shape" scaling 
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functions which are not in accordance with data analy- 
sis [35, 44, 47]. 

In this work wc extend the analysis of scaling to the 
relativistic GF approach, whose predictions will be com- 
pared with previous results obtained within the RMF, 
RPWIA, and rROP frameworks. This study will show to 
which degree superscaling is fulfilled by the GF calcula- 
tions, and, moreover, how the GF scaling function com- 
pares with the experiment. As already mentioned, the 
GF treats FSI consistently in the inclusive and exclusive 
reactions, hence its application to scaling and superscal- 
ing ideas, which emerge as an essential outcome of nature, 
results mandatory. This allows us to get a clear insight 
concerning the capability of the relativistic GF approach 
applied to QE inclusive (e, e') reactions, as well as the 
uncertainties ascribed to the various ingredients of the 
model, particularly the specific energy-dependent terms 
in the complex optical potentials involved. 
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FIG. 1: Differential cross section of the ^^C{e, e') reaction for 
an incident electron energy e = 1 GeV and three values of the 
momentum transfer, i.e., q = 500, 800, and 1000 McV/c, cal- 
culated by the Pavia (solid) and the Madrid-Sevilla (dashed) 
groups with RPWIA (left) and rROP (right). 



III. RESULTS AND DISCUSSION 

In this section the numerical results of the different rel- 
ativistic models developed by the Pavia and the Madrid- 
Sevilla groups to describe FSI in the inclusive quasielastic 
electron-nucleus scattering are considered. As a first step 
results obtained by the two groups in the RPWIA and 
rROP approaches are compared in order to check the 
consistency of the numerical programs when calculations 
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FIG. 2: Differential cross section of the ^^C(e, e') reaction 
in the same kinematics as in Fig. 1. The solid and long 
dot-dashed lines are the GF results calculated with the two 
different optical potentials EDADl and EDAD2. The dashed 
lines are the results of the RMF model. 



are carried out under the same conditions. Then the 
results corresponding to the RMF model developed by 
the Madrid-Sevilla group and the relativistic GF model 
developed by the Pavia group are compared. This com- 
parison is performed for the ^^C(e, e') cross section and 
scaling function calculated with various models under dif- 
ferent kinematics. 



A. Differential cross section 

The ^^C (e, e') cross section has been calculated for a 
fixed value of the incident electron energy e = 1 GeV 
and three values of the momentum transfer, i.e. q = 
500, 800 and 1000 MeV/c. The relativistic initial states 
are taken as Dirac-Hartree solutions of a relativistic La- 
grangian written in the context of a relativistic mean field 
theory with the NLSH parameterization [79-81]. 

The cross sections calculated in RPWIA and in rROP 
by the Pavia and the Madrid-Sevilla groups are com- 
pared in Fig. 1. The calculations are performed using 
the same ingredients for the relativistic wave-functions 
and the current operator. Almost identical results are 
obtained in RPWIA. In rROP the two results are very 
similar, up to a few percent. 

The comparison in Fig. 1 is a first important and nec- 
essary benchmark of the two independent computer pro- 
grams, which allows us to estimate the numerical uncer- 
tainties and gives enough confidence on the reliability of 



6 



q = 500 MeV/c 




GF1 






GF2 _ 

RMF 






rROP 


ii 1 <■': 
- i f': 
\' :' 


\^ 


RPWIA 


1 .• ^ 







300 400 

0) [MeV] 




(J [MeV] 

FIG. 3: Differential cross section of the ^^C(e, e') reaction 
for e = 1 GeV and q = 500 and 1000 MeV/c. The sohd, 
long dot-dashed, and dashed lines are the same as in Fig. 1. 
The dot-dashed and dotted lines are the rROP and RPWIA 
results, respectively. 
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FIG. 4: Differential cross section of the ^^C(e, e') reaction for 
different beam energies and electron scattering angles. Line 
convention as in Fig. 2, experimental data from [85-87]. 



both calculations. Further discrepancies between the re- 
sults of the two groups can now be ascribed to differences 
of the models but not to inconsistencies in the calcula- 
tions. 

The cross sections evaluated with the RMF and GF 

models for the same kinematics as in Fig. 1 are presented 
in Fig. 2. In the case of the GF approach two differ- 
ent parameterizations for the rclativistic optical potential 
have been used: the energy-dependent and A-dependent 
EDADl and EDAD2 complex phenomenological poten- 
tials of [82], which are fitted to proton elastic scattering 
data on several nuclei in an energy range up to 1040 
MeV. In the figures the results obtained with EDADl 
and EDAD2 arc denoted as GFl and GF2, respectively. 
The differences between the results of the RMF and GF 
models increase with the momentum transfer. Also dis- 
crepancies between the GFl and GF2 cross sections de- 
pend on the momentum transfer. kX. q = 500 MeV/c the 
three results in Fig. 2 are similar, both in magnitude and 
shape. At larger q, q = 800 MeV/c, moderate differences 
are found, whereas the discrepancy between the three 
approaches gets larger at g = 1000 MeV/c. The shape 
of the RMF cross section shows an asymmetry, with a 
long tail extending towards higlica- values of w, that is 
essentially due to the strong energy-independent scalar 
and vector potentials present in the RMF approach. On 
the contrary, in the case of GFl and GF2, the asymme- 
try towards higher w is less significant but still visible. 
The GFl and GF2 cross sections show a similar shape 



but with a significant difference in the magnitude. At q 
= 1000 MeV/c both of them are higher than the RMF 
cross section in the region where the maximum occurs. 
However, note that a stronger enhancement is obtained 
with GFl, which at the peak overshoots the RMF cross 
section up to 40%. Overall, taking into account the dif- 
ferent ingredients that these estimates of the inclusive 
response contain, the similarity of the predictions for the 
inclusive cross-section is remarkable, particularly at the 
two lower values of q. The larger differences seen for 
the largest q value, not only between the RMF and GF 
models, but also between the two GF results, is simply 
an indication of the difference in the ingredients of these 
calculations. 

Indeed, the different behaviour presented by the RMF 
and GF results as a function of q and uj is directly linked 
to the specific structure of the relativistic potentials in- 
volved in the RMF and GF models. Whereas the RMF is 
based on the use of a strong energy-independent real po- 
tential, the GF approach makes use of a complex energy- 
dependent optical potential. In GF calculations the be- 
havior of the optical potential changes with the momen- 
tum and energy transferred in the process, and higher 
values of q and lo correspond to higher energies for the 
optical potential. The results obtained with GFl and 
GF2 are consistent with the general behaviour of the 
phenomenological rclativistic optical potentials and are 
basically due to their imaginary part. To make this clear, 
we present again in Fig. 3 the GF and RMF cross sec- 
tions for q = 500 and 1000 MeV/c compared directly 
with the corresponding results obtained within the RP- 
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FIG. 5: Longitudinal contribution to the scaling function for 
three values of the momentum transfer, i.e., q = 500, 800, and 

1000 McV/c, obtained by the Pavia (solid) and the Madrid- 
Sevilla (dashed) groups with RPWIA (left) and rROP (right). 
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FIG. 6: Longitudinal contribution to the scaling function for 
q = 500, 800, and 1000 MeV/c with the GFl (solid), GF2 
(long dot-dashed), and RMF (dashed) models compared with 
the averaged experimental scaling function. 



WIA and rROP models. It is known that the real terms 
of the relativistic optical potentials are very similar for 
the different parameterizations. This explains why the 
cross section evaluated within the rROP approach docs 
not show sensitivity to the particular parameterization 
considered for the ROP. On the other hand, the energy- 
dependent scalar and vector components of the real part 
of the ROP get smaller with increasing energies. Thus 
the rROP result approaches the RPWIA one for large 
values of cj. In contrast, the imaginary (scalar and vec- 
tor) part presents its maximum strength around 500 MeV 
being also sensitive to the particular ROP parameteriza- 
tions. This explains the differences observed between the 
rROP and the two GF results as a function of ui and q. 
This significant discrepancy between GF and rROP cross 
sections in Fig. 3 seems to contradict previous results 
shown in [37]. However, kinematical conditions reported 
in [37] corresponded to lower values of the momentum 
transfer, where there is no reason a priori to expect the 
rROP and GF predictions to be closer or further apart. 

Of particular relevance is the difference between the 
GFl and GF2 results. These are obtained with optical 
potentials that reproduce the elastic proton-nucleus phe- 
nomenology to a similar degree [82]. However, one must 
keep in mind that elastic observables do not completely 
constrain optical potentials, and indeed it has been of- 
ten seen how the predictions of the EDADl and EDAD2 
potentials for nonelastic observables, such as for instance 
(e, e'p) or electron trasparencies [56, 57], differ signifi- 
cantly. The differences between GFl and GF2 are mostly 



due to their different imaginary part, that includes the 
overall effect of the inelastic channels and is not uni- 
vocally determined by the elastic phenomenology. The 
most convenient choice of the phcnomelogical optical po- 
tential to be employed within the GF approach should 
thus be made after a comparison to inclusive data. 

In Fig. 4 the GFl, GF2, and RMF results are compared 
with the experimental cross sections for three different 
kinematics [85-87]. A recent review of the experimen- 
tal situation as well as different theoretical approaches 
can be found in [2]. Results in Fig. 4 show that the 
three models lead to similar cross sections. The main 
differences are presented for higher values of the momen- 
tum transfer, about 800 MeV/c (bottom panel), where 
the GFl cross section (solid line) is larger than the GF2 
(dot-dashed) and RMF (dashed) ones. The experimental 
cross section is slightly underpredicted in the top panel 
and well described in the middle panel by all calculations. 
Finally, results in the bottom panel show a fair agreement 
with data in the case of GFl, whereas GF2 and RMF 
underpredict the experiment. Summarizing, the compar- 
ison with data, although satisfactory on general grounds, 
gives only an indication and cannot be conclusive until 
contributions beyond the QE peak, like meson exchange 
currents and Delta effects, which may play a significant 
role in the analysis of data even at the maximum of the 
QE peak, are carefully evaluated [20, 45, 46]. These pro- 
cesses contribute to the cross-section and the comparison 
of the pure nuclconic predictions of the GFl, GF2, and 
RMF models to data will only indicate what is the con- 
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FIG. 7: Longitudinal contributions to the scaling function 
for q = 500 and 1000 MeV/c compared with the averaged 
experimental scaling function. Line convention as in Fig. 3. 
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FIG. 8: Analysis of first-kind scaUng, fLii^') for q — 500 
(solid), q = 800 (dot-dashed), and 1000 MeV/c (dashed) with 
the GFl (top panel), GF2 (middle panel), and RMF (bottom 
panel) models using the same results already displayed in Fig. 
6. 



tribution of the non-nucleonic degrees of freedom to the 
cross-sections. 



B. Scaling functions 

The effects already discussed for the differential cross 

sections are also present in the scaling functions. Here 
we compare results for the longitudinal component of the 
scaling function JlH^') using the same descriptions for 
the final state interactions already considered for the dif- 
ferential cross sections. 

As a first step, the longitudinal contribution /l(V'') 
obtained in RPWIA and rROP by the Pavia and the 
Madrid-Sevilla groups at three values of the momentum 
transfer is displayed in Fig. 5 showing almost coincident 
results. Similar comments apply to the transverse contri- 
bution /t(V'')- These results, in addition to the cross sec- 
tions in Fig. 1, confirm the consistency of the numerical 
codes when calculations are performed under the same 
conditions. 

The scaling function JlHj') evaluated within RPWIA 
and rROP shows a very mild dependence on the momen- 
tum transfer for both positive and negative values of the 
scaling variable ip', i.e., violation of scaling of first kind 
is small. 

In Figs. 6 and 7 we compare evaluated with 

different models and for several values of the momentum 
transfer with the averaged QE phenomenological scaling 

function extracted from the analysis of (e, e') data [41- 
43]. As already shown in previous works [35, 44, 47], the 



RMF model produces an asymmetric shape with a long 
tail in the region with tp' > that follows closely the 
phenomenological function behavior. The asymmetry in 
the data has usually been attributed to physical effects 
beyond the mean field such as two-body currents and 
short-range correlations [24]. Within a non-relativistic 
framework such contributions need to be considered in 
order to get asymmetry [24, 88, 89]. By contrast, the 
RMF approach is capable of explaining the asymmetric 
behavior of data within the framework of the relativis- 
tic impulse approximation taking advantage of its strong 
relativistic scalar and vector potentials. The results with 
the GF model are similar to those obtained with RMF 
a.t q = 500 MeV/c and, with moderate differences, at q 
— 800 MeV/c, while visible discrepancies appear at q = 
1000 MeV/c. On the other hand, discussion of results 
for the scaling functions follows similar trends to the one 
already applied to the behaviour of the cross sections in 
Fig. 2, i.e., at higher q- values the maximum strength oc- 
curs for the GFl model being the RMF one the weakest. 

The asymmetric shape with a tail in the region of pos- 
itive ip' is obtained in both RMF and GF models which 
involve descriptions of FSI either with a strong energy- 
independent real potential or with a complex energy- 
dependent optical potential, respectively. The scaling 
functions corresponding to RPWIA and rROP, which are 
also presented in Fig. 7 do not show any significant asym- 
metric tail for ip' > 0. The different dependence on the 
momentum transfer shown by the potentials involved in 
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the RMF and GF approaches makes the GF scaling func- 
tion tail less pronounced as the value of q goes up. 

The comparison of the different models to the longi- 
tudinal scaling function is illuminating. Wc must recall 
that the experimental longitudinal response can be con- 
sidered as a much better representation of the pure nu- 
cleonic contribution to the inclusive cross-section than 
the total cross-section. It is remarkable that, as seen in 
figs. 6 and 7, except for the highest value of q considered 
(1000 MeV/c), GFl, GF2 and RMF approaches yield 
very similar predictions for the longitudinal response, in 
good agreement with the experimental longitudinal re- 
sponse. The asymmetric tail of the data and the strength 
at the peak are fairly reproduced by the three approaches. 
However, for q = 1000 MeV/c, only the RMF approach 
seems to be favoured from the comparison to the data, 
while GFl and GF2 yield now rather different predic- 
tions than the RMF approach, that seem to be ruled out 
by data. We must keep in mind that the GF approach 
uses as input the phenomenological optical potentials. It 
is clear that, as the momentum transfer increases, the 
phenomenological optical potential will (implicitely) in- 
corporate a larger amount of contributions from non nu- 
cleonic degrees of freedom, such as for instance the loss 
of (elastic) flux into the inelastic delta excitation with or 
without real pion production. Thus the input of the GF 
formalism is contaminated by non purely nucleonic con- 
tributions. Consequently, GF predictions depart from 
the experimental QE longitudinal response, that effec- 
tively isolates only nucleonic contributions. This differ- 
ence, which is larger with increasing momentum transfer, 
emerges as an excess of strength predicted by the GF 
model as it translates a loss of flux due to non-nucleonic 
processes into inclusive purely nucleonic strength. On 
the other hand, the RMF model uses as inpiit the effec- 
tive mean fleld that reproduces saturation properties of 
nuclear matter and of the ground state of the nuclei in- 
volved, and thus it is much more suited to estimate purely 
nucleonic contribution to the inclusive cross-section, even 
at q = 1000 MeV/c. Taking these facts into account, 
the comparison of the models to the data, both for total 
cross-sections and longitudinal responses yields what one 
reasonably expects. 

An analysis of the scaling of flrst-kind, i.e., indepen- 
dence of the momentum transfer, is illustrated in Fig. 8. 
The results are the same already shown in Fig. 6, but are 
presented in a different way. Each panel corresponds to 
a specific description of FSI (GFl, GF2, and RMF) and 
includes results obtained for different values of q. The 
experimental data are compatible with a mild violation 
of the first-kind scaling, particularly in the positive if)'- 
region. In Refs. [35, 44] the scaling functions evaluated 
with the RPWIA and rROP models were shown to de- 
pend very mildly on the transferred momentum in the 
whole, positive and negative, ip' region. In the case of 
the RMF approach, there is a slight shift in the region 
■0' < 0, whereas the model breaks scaling approximately 
at 30% level when > Q. Similar results are obtained 



with the GF models, where a shift in the region of nega- 
tive ij)' also occurs, and scaling is broken for tl)' > 0. This 
scaling violation for ijj' > is larger with GFl due to the 
enhancement produced in the region of the QE peak by 
this phenomenological optical potential at higher values 
of q. As a consequence, the comparison between the ex- 
perimental QE scaling function and the GFl results is 
worse than the corresponding comparison for GF2. 



IV. SUMMARY AND CONCLUSIONS 

The work developed in this paper has emerged as a 
close collaboration between the Pavia and Madrid-Sevilla 
groups. For the last years, both groups have been deeply 
involved in studies of lepton scattering reactions with nu- 
clei. Inclusive and exclusive electron scattering processes 
have been analyzed, as well as reactions involving the 
use of neutrinos. The calculations performed by the two 
groups, which are similar in some aspects, i.e., treatment 
of relativistic ingredients, bound-nucleon wave functions, 
single- nucleon current operators, etc., present also clear 
differences concerning the description of final state in- 
teractions, which constitute an essential ingredient for a 
successful comparison with data. 

The consistency of the numerical calculations devel- 
oped by our two groups is checked by comparing results 
in the plane wave limit (RPWIA) and making use of the 
real part of the relativistic optical potential (rROP). Al- 
most identical predictions come out within RPWIA and 
very similar ones with rROP. This reinforces our confi- 
dence on the reliability of both calculations. As known, 
the description of inclusive (e,e') reactions requires the 
contribution from the inelastic channels to be retained. 
Within the framework of the relativistic impulse approxi- 
mation, a simple recipe to compute the inclusive strength 
is the use of purely real potentials in the final state. This 
is the case of the rROP approach (phenomenological rela- 
tivistic optical potential but without the imaginary part). 
However, although rROP conserves the flux, its use for 
inclusive processes is not entirely satisfactory since the 
(complex) relativistic optical potential has its origin in 
the description of exclusive processes where only the clas- 
tic channel contributes to the observables. Two other 
models based on the relativistic impulse approximation 
have been used recently to account for FSI. The former 
employes distorted waves obtained with the same rela- 
tivistic mean field considered for the initial bound states. 
This is denoted as RMF and it fulfills the dispersion rela- 
tion and the continuity equation. The latter procedure is 
based on the relativistic Green's fimction (GF) approach, 
which allows one to treat FSI consistently in the inclusive 
and exclusive reactions. 

Differential cross sections and scaling functions eval- 
uated with both models for different kinematical situa- 
tions have been compared. Discrepancies are shown to 
increase with the momentum transfer. This is linked 
to the energy-dependent optical potentials involved in 
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the relativistic GF method by contrast to the strong 
energy-independent RMF potentials. Moreover, results 
presented for two different parameterizations of the ROP, 
namely, EDADl vs EDAD2, prove the importance of 
the imaginary term, which gets its maximum strength 
around 500 MeV, whereas the real part gets smaller as 
the energy increases. This explains the different behavior 
shown by the results of the GF and rROP models, the 
latter approaching the RPWIA results for higher values 
of the transferred energy and momentum. Furthermore, 
discrepancies between the GF and RMF results are also 
clearly visible as q goes up, due to the contribution from 
non-purely nucleonic inelasticities to the phenomenolog- 
ical optical potential 

All models considered respect scaling and superscal- 
ing properties. Furthermore, the significant asymmetry 
in the scaling function produced by the RMF is strongly 
supported by data [35]. On the contrary, asymmetry is 
largely absent in RPWIA and rROP predictions. The 
relativistic GF approach leads to similar results to RMF, 
i.e., with presence of the asymmetry for intermediate q- 
values. On the contrary, visible discrepancies emerge 
for larger q, being the GF scaling fimction tail less pro- 
nounced but showing more strength in the region where 
the maximum occurs. Moreover, the GF results for high 
q present a strong dependence on the specific parameteri- 
zation considered for the optical potential. These results 
are directly linked to effects introduced by the imagi- 
nary (scalar and vector) term in the optical potential that 
presents a high sensitivity to the particular ROP param- 
eterization. The relativistic GF approach, even based 
on the use of a complex optical potential, is shown to 
preserve flux conservation, hence the imaginary term in 
the potential leading to a redistribution of the strength 
among different channels. This explains the difference 
observed between RMF and GF predictions, the latter 
with additional strength in the region close to the max- 
imum in the QE response. This behavior could be con- 
nected with effects coming from the contribution of the 
A which are, somehow, accounted for in a phenomeno- 
logical way by the GF approach, modifying consequently 



the responses even in the region where the QE peak gives 
the main contribution. Notice that the higher the trans- 
fer momentum is, a stronger overlap between the QE and 
A peaks occurs. This makes very difficult to isolate con- 
tributions coming from either region. 

Although great caution should be exercised in extend- 
ing the above comments before more conclusive studies 
are performed, the present analysis can be helpful in dis- 
entangling different treatments of FSI and their connec- 
tions with different physics aspects involved in the pro- 
cess. The similarities of the GF and RMF predictions for 
the inclusive cross-sections, particularly for intermediate 
values of g, in spite of the very different phenomenological 
ingredients they consider, and the very reasonable agree- 
ment with the data for the longitudinal scaled response, 
that constitutes a good representation of the experimen- 
tally measured purely nucleonic response to the inclusive 
cross-sections, are a clear indication of the fact that both 
models make a very decent job in estimating the inclu- 
sive contribution. It will be interesting to investigate 
the possibility of disentangling in the phenomenological 
optical potential the contributions due to non-nucleonic 
inelasticities and extract a 'purely nucleonic' optical po- 
tential which could then be used in the GF approach and 
contrasted against the experimental longitudinal scaling 
function. This work can be considered as a first step in 
this direction. 
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